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ABSTRACT 

The results of a detailed investigation of the distribution of stresses in aluminum and 
composite panels subjected to uniform end shortening are presented. The focus problem is a 
rectangular panel with two longitudinal stiffeners, and an inner stiffener discontinuous at a 
central hole in the panel. The influence of the stiffeners on the stresses is evaluated 
through a two-dimensional global finite element analysis in the absence or presence of the 
hole. Contrary to the physical feel, it is found that the maximum stresses from the global 
analysis for both stiffened aluminum and composite panels are greater than the corresponding 
stresses for the unstiffened panels. The inner discontinuous stiffener causes a greater 
increase in stresses than the reduction provided by the two outer stiffeners. A detailed 
layer-by-layer study of stresses around the hole is also presented for both unstiffened and 
stiffened composite panels. A parallel equation solver is used for the global system of 
equations since the computational time is far less than that using a sequential scheme. A 
parallel Choleski method with up to 16 processors is used on Flex/32 Multicomputer at NASA 
Langley Research Center. The parallel computing results are summarized and include the 
computational times, speedups, bandwidths, and their inter-relationships for the panel 
problems. It is found that the computational time for Choleski method decreases with a 
decrease in bandwidth, and better speedups result as the bandwidth increases. 

1 Professor, Department of Civil Engineering, Old Dominion University, Norfolk, Virginia 
23529. Fellow, ASCE. 

2 Graduate Student, Department of Civil Engineering, Old Dominion University, Norfolk, 
Virginia 23529. 
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I. INTRODUCTION 


There is an ever increasing need for accurate and high speed numerical computations to 
evaluate the stress distributions in aerospace structures. Furthermore, the constant demand 
to reduce the weight/stiffness ratios for such structures has made the use of composite 
materials an attractive option in comparison to aluminum alloys. The process of routine 
stress analysis as well as that of reducing the structural weight is considerably improved 
through a detailed global and/or local stress analysis of structural components in the 
presence or absence of discontinuities. This report presents the outcome of a two- 
dimensional sequential and concurrent finite element stress analysis of unstiffened and 
stiffened rectangular panels with and without a central hole. Both aluminum and composite 
panels are investigated. Such panels have wide applications in wing and fuselage structures 
of aircraft. Holes are drilled into panels used in aerospace vehicles to facilitate bolting and 
riveting to primary structure, or to provide access to the interior of wings. 

The buckling and failure mechanisms for rectangular panels with and without stiffeners 
and a central hole are studied by some researchers. Examples include the experiments by 
Starnes and Williams (reference 1), and Knauss and Henneke (reference 2) on composite 
panels. A limited number of three-dimensional finite element analyses of unstiffened 
rectangular panels with hole are available. References 3 through 7 are representative 
examples published in over a decade, and contain results for composite panels subjected to 
uniform compressive or tensile boundary stresses. Hirai, Wang, and Pilkey (reference 8) 
presented a finite element zooming technique employing static condensation and structural 
reanalysis. They have reported stress concentration factors with various types of local 
meshes for a rectangular steel panel with hole and subjected to a uniaxial tensile boundary 
stress. A review of the existing literature indicates that no detailed pre-buckling numerical 
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stress distributions are published in the past for rectangular panels with stiffeners and a 
central hole. 

The objective of this report is to study and compare the distribution of stresses and 
maximum stresses near the hole in unstiffened and stiffened, aluminum and composite panels, 
without and with a hole, when subjected to uniform end shortening. The maximum stress 
values based on global analysis for both stiffened aluminum and composite panels are found 
to be greater than the corresponding values for unstiffened ones. This phenomenon is 
contrary to the common belief that the stiffeners reduce such stresses. Hence, a further 
investigation became necessary, and was accomplished through the identification of the 
individual effects of the inner and outer stiffeners on the maximum stresses. Also, only a 
limited amount of work has been carried out in the past on finite element stress analysis 
problems through the use of computationally efficient parallel solvers for simultaneous 
equations. The performance of a parallel equation solver based on Choleski method 
(references 9 and 10) is also evaluated for the panel problems considered herein. 
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II. DESCRIPTION OF PROBLEMS 


Figure 1 shows six types of panels analyzed in this paper. Figures 1(a) and 1(b) show 
rectangular panels without, and with stiffeners, respectively, in the absence of a hole. The 
panels shown in Figures 1(c) through 1(f) have a central hole. The panel in Figure 1(c) is 
unstiffened, while those in Figures 1(d) through 1(f) are provided with one, two, and three 
stiffeners, respectively, as shown. The problems shown were solved using the NASA Langley 
Computational Structural Mechanics testbed software system (reference 11). 

The finite element discretizations for the panel skin and stiffeners are shown in Figure 
2. The overall panel dimensions Lj and L 2 ; the distance between panel edge and outer 
stiffener, B e , and the distance between the inner and outer stiffeners, B s , are shown in 
Figure 2(a). Each stiffener has a height of H s as indicated in Figure 2(b). In these figures, 
the global axes are indicated as Xj, X 2 , and X 3 . The discretization in the region around the 
hole used in the global analysis is shown in Figure 2(c). The global model in the region 
around the hole, and the local model around the hole are shown in Figure 2(d). The 
discretization for the stiffened panel with a hole is shown in Figure 2(e). Quadrilateral 
elements are used for the panels considered, except for the panels without a hole in which 
case additional triangular elements are used to fill in the region of the hole. Each node has 
two translational degrees-of-freedom, along the Xj and X 2 axes, for the unstiffened panels. 
However, for the stiffened panels, each node has six degrees-of-freedom, namely, three 
translations and three rotations relative to the Xj, X 2 , and X 3 axes. 

The following dimensions are adopted for the panel skin and stiffeners: Lj = 30.0 in.; L 2 
= 11.5 in.; R = 1.0 in.; B s = 4.5 in.; B e = 1.25 in.; H s = 1.4 in., and A = 4.0 in. The thickness 
of the aluminum panel is taken as 0.1 in. for both the panel skin and each stiffener. The 
thickness of the composite panel skin is 0.1375 in. for a total of 25 plies, and that of each 
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composite stiffener is 0.132 in. for a total of 24 plies. The Young’s modulus, and Poisson’s 
ratio for aluminum are taken as 10,000 ksi, and 0.30, respectively. The material of the 
composite panels is graphite-epoxy (T300/5208) with the following mechanical properties: Ej = 
19,000 ksi; E2 = 1,890 ksi; = 0-38, and Gj 2 = 930 ksi, as adopted in reference 2. 

Each panel is subjected to a uniform edge shortening of 0.1 in. in the longitudinal (Xj) 
direction, imposed at the face OS as shown in Figure 2(a). The boundary conditions adopted 
for the unstiffened, and stiffened panels are summarized in Tables 1, and 2, respectively. 
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III. STRESS DISTRIBUTIONS 


In this section, the results are given for the various types of panels shown in Figure 1. 

Dimensionless contour plots for the principal compressive stresses are presented for a quarter 

of these panels, based on the global analysis. Contours are also drawn for the region around 

the hole for the unstiffened panel, based on the local analysis. For the stiffened panel, the 

maximum midsurface shear stress distribution at the location of the stiffeners is also given. 
* 

Both aluminum and composite panels are considered. Furthermore, the influence of inner and 
outer stiffeners on the stresses is identified. 

3.1 Stress Distributions in Aluminum Panels 

3.1.1 Unstiffened and Stiffened Aluminum Panels without Hole 

The total degrees-of-freedom (DOF) for unstiffened, and stiffened panels without hole 
are 554, and 2346, respectively. Table 3 presents the dimensionless principal stresses, 02 , 
given by the ratio of the principal compressive stress, a 3 , and the far field (boundary) 
stress, a 0 , for the panels shown in Figures 1(a) and 1(b). The values are summarized in the 
table for various locations identified by the coordinates Xj, X 2 indicated in Figure 2. The 
corresponding stresses in the stiffeners are summarized in Table 4. 

3.1.2 Unstiffened and Stiffened Aluminum Panels with Hole 

The 02 contours for the unstiffened panel with hole based on a 552 DOF global model 
are presented in Figure 3(a). The contours based on the local analysis for the region around 
the hole are shown in Figure 3(b). The local analysis procedure is based on imposing an 
interpolated displacement field for the region around the hole . The stress concentration 
factors (SCF) for the global and local analysis are 2.74, and 3.00, respectively. Figure 4(a) 
shows the stress contours for the stiffened panel skin analyzed globally in which the outer 
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and inner stiffener locations are also indicated. Figures 4(b) and 4(c) show the contours for 
these stiffeners. The resulting SCF for the stiffened aluminum panel is 2.89. 

It is interesting to note that the SCF for the stiffened panel is gre.ater than that of 
the unstiffened one. The reason for this phenomenon is explained later in this report. The 
outer stiffeners develop higher stresses than those in the inner one. Also, a substantial 
part of the inner stiffener is stress-free. The midsurface maximum shear stress (»- max ) 
distribution for the panel at the location of both inner and outer stiffeners is shown in 
Figure 5 as a function of the distance from the center of the hole. The r max exhibits a 
greater variation along the panel length at the inner stiffener location in comparison to that 
observed for the outer stiffener location. The reason may be attributable to a more 
dramatic variation in normal stresses for the inner stiffener in the panel longitudinal 
direction as compared to that for the outer stiffeners, as may be seen readily from the 
nature of the contours in Figures 4(b) and 4(c). 

3.2 Stress Distributions in Composite Panels 

The results for various types of composite panels shown in Figure 1 are presented in 
this section. Dimensionless contours for the principal compressive stress resultants based on 
global analysis are presented for a quarter of each panel shown in Figures 1(c) and 1(f). 
Contours are also drawn for the region around the hole for the unstiffened panel based on 
the local analysis. Tables showing the variation of principal stresses and maximum shear 
stresses around the hole and through the 25 layers are presented for both stiffened and 
unstiffened panels. Also included are the plots showing the variation of principal stresses 
and maximum shear stresses around the hole in the layer with highest value of stresses. In 
addition, plots are presented for the variation of these stresses through the layers and at the 
node with maximum stresses. 


6 



3.2.1 Unstiffened and Stiffened Panels without Hole 

Table 5 presents the dimensionless principal stress resultant, 5 given by the ratio of 
the principal compressive stress resultant, c^, and the far field (boundary) stress resultant, 
<r sro> for the panels shown in Figures 1(a) and 1(b). The stress resultant, c^, is defined by: 

o = At Efj 

sr 1 

in which At is the thickness of one ply, fj is the stress in the i-th ply, and the summation 
is over N plies. The dimensionless stress resultants are summarized in the table for various 
locations identified by the coordinates X\ and X 2 indicated in Figure 2. The corresponding 
stress resultants in the stiffeners are summarized in Table 6. It is seen that in the 
unstiffened panel a : . ranges from 1.00 to 1.13, while in the stiffened panel, it ranges from 
0.96 to 1.05. The corresponding values for the stiffener are in the range from 1.09 to 1.39. 

3.2.2 Unstiffened Panel with Hole 

Figure 6(a) shows the dimensionless stress resultant o ^ contours for the 552 DOF 
global model. The contours based on the local analysis for the region around the hole are 
shown in Figure 6(b). The maximum 5^ value for the global and local analysis are 1.74 and 
1.90, respectively. Table 7 presents the variation of maximum principal stress, <x\, around the 
hole for all the 25 layers of the panel skin. It is seen that the a\ value is maximum at d" 
for a -45° ply-orientation. Figures 9 and 10 show the variation of a\ around the hole, and 

It H 

through the layers at D , respectively. The maximum value of <tj exists at D , while the 
minimum values are at C and c". It is also seen that the same values of a\ exist for a 
given orientation through the layers at each node. Table 8 presents the variation of 
minimum principal stress, 02 , around the hole for all the 25 layers of the panel skin. The 
maximum value is seen at C for 0’ ply-orientation. Figure 11 and 12 show the variation of 


7 



ai around the hole, and through the layers at C\ respectively. It is seen that the maximum 

values of <*2 are at C and C , and the same values of *2 exist for a given orientation 

through the layers at each node. Table 9 shows the variation of maximum shear stress, 

r , around the hole for all the 25 layers of the panel skin. Here too, the maximum value 
max’ 3 

is at C for 0° ply-orientation. Figures 13 and 14 show the variation of t around the 

max 

hole, and through the layers, respectively. The maximum values are at C and c" and the 

same values of r exist for a particular ply-orientation through the layers at each node, 
max 

3.2.3 Stiffened Panel with Hole 

The dimensionless stress resultant 5^ contours for the panel skin of the 2328 DOF 

global model are presented in Figure 7(a). The stress resultant contours for the stiffeners 

are shown in Figures 7(b) and 7(c). The maximum 5^ f° r the global analysis is found to be 

1.8. The distribution of the maximum shear stress resultants r , for the inner and 

max 

outer stiffeners is shown in Figure 8. A rapid variation is found at the inner stiffener 
location except for the region close to the panel edge, that is, as approaches 15 in. Its 
variation at the outer stiffener location is relatively small. Table 10 shows the variation of 
maximum principal stress, o\, around the hole for all the 25 layers of the panel skin. The 
maximum value of <tj is seen at d" on the bottom-most layer (25th) with a ply-orientation of 
45°. Figure 15 shows the variation of o\ around the hole for the 25th layer, and Figure 16, 
the variation through all the layers at d". The highest value of is at D , while the 

i* 

lowest value is at C . Table 11 presents the variation of minimum principal stress, *2, 

around the hole for all the 25 layers of the panel skin. The maximum value of 02 is seen at 

C' on the third layer, with a ply-orientation of 0°. Figure 17 shows the 02 variation around 

the hole for the 3rd layer, and Figure 18, the variation through all the layers at C'. The 

maximum value of 02 is at C\ while the minimum is at d". Table 12 shows the variation of 

maximum shear stress, r around the hole for all the 25 layers of the panel skin. The 

max, 
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maximum value of r is found at d" on the 25th layer with a 45° ply-orientation. Figure 
max 

19 shows the r variation around the hole for the 25th layer, and Figure 20 the variation 
max 

through all the layers at D . The maximum value of r is at D , while the minimum is at 

max 

a location half-way between c" and D*. 

3.2.4 Nature of Stresses in Composite Panels 

The unstiffened panel without a hole exhibits a fairly uniform stress resultant 

throughout, except at the center of the panel where it is nearly 13% more than the 

surrounding stress resultants. In the stiffened composite panel, a nearly uniform stress 

resultant distribution in the range from 0.95 to 1.05 of the boundary stress resultant is seen. 

The stiffeners exhibit larger values of upto 1.39. In the unstiffened panel high stress 

gradients exist in the vicinity of the hole and are more prominent along C'C as seen in the 

Figure 7. The stiffened panel with a hole has higher stress resultants at C’ than that in the 

unstiffened one. The outer stiffeners develop higher stress resultants than the inner 

stiffeners. Furthermore, as observed in the case of aluminum stiffener, the composite 

stiffener also possesses a region with little or no stress resultant. The r variation is 

max 

more in the inner stiffener location than at the outer stiffener location due to considerable 

variation of normal stress resultants in the panel longitudinal direction. 

The maximum values of principal stresses a[ and <?2 occur on different layers for both 

unstiffened and stiffened composite panels with a hole. The unstiffened panel has the 

highest values of ^ and maximum shear stress, r , on the same layer, while the stiffened 

max 

panel has the maximum values of at and r on the same layer. The stresses are the same 

1 max 3 

for a particular ply-orientation for the unstiffened panel whereas they are different for each 
layer for the stiffened one. 

For both the unstiffened and stiffened panels, the a\ value is maximum at the location 

n 

D , but it exists on the -45° ply-orientation for the former and 45° for the latter. The a2 
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distribution around the hole for the unstiffened panel shows a symmetry about d"d i . The 

maximum value of ^2 * s at C" and c" while its minimum value is at d" for both the cases. 

For the stiffened panel, the a 2 distribution through the layers has its maxiipum values on the 

third layer, with a 0° ply-orientation. The same is true for the unstiffened panel except 

that it is repeated for each ply with 0° orientation. The r distribution around the hole 

max 

n I f 

is symmetrical about D D for only the unstiffened panel with its maximum value at C and 
C on the third layer, and repeated for each layer with the same orientation. The stiffened 
panel develops the maximum value in the 25th layer at d". 


3.3 Influence of Inner and Outer Stiffeners on Stresses 

As mentioned previously in this report, the SCF and ( 5 sr 2) max va l ues for stiffened 

aluminum and composite panels, respectively, were found to be greater than the 

corresponding values for the unstiffened ones. To understand the cause of this phenomenon 

global analyses of the panels shown in Figures 1(d), and 1(f) were conducted while utilizing 

1,920 DOF and 2,070 DOF model, respectively. The resulting SCF and (5 _) values for 

srz max 

aluminum and composite panels, respectively, together with those obtained previously for the 
panels shown in Figures 1(a) and 1(f), are summarized in Tables 13 and 14. 

The results in Table 13 reveal that when an aluminum panel is "stiffened* by the inner 
stiffener only (Case 2), the SCF is 2.93, that is about 7% greater than that of the 
unstiffened panel (Case 1). The inner stiffener appears to be acting as a nearly direct load- 
transfer unit between the loaded boundaries and the hole thereby causing a further increase 
in the maximum compressive stresses. The SCF value for the aluminum with all three 
stiffeners (Case 4) is found to be 2.89 which is also greater than that for the unstiffened 
one. However, a comparison of the SCF values for Cases 1,2 and 4 leads to the conclusion 
that the outer two stiffeners must be introducing a beneficial effect on the maximum 
stresses. This conclusion was confirmed through the analysis of the panel with only the 
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outer stiffeners (Case 3) which provided an SCF of 2.69, that is, lower than that for the 
unstiffened panel. 

The results in the Table 14 reveal that for the panel with only the middle stiffener (case 
2), the maximum dimensionless stress resultant, (3 -) , is 1.84 or 6% greater than that 

for the unstiffened panel. This increase is attributed to the middle stiffener acting as a 
direct load-transfer unit between the loaded boundaries and the hole causing a further 
increase in the maximum compressive stress resultants. The beneficial effect of introducing 


two outer stiffeners (Case 3) is evident from 


(7 -) value of 1.66, which is 5% less than 
v sr2 / max 


that for the unstiffened one. The adverse effect of the middle stiffener is also seen for the 


panel with all the three stiffeners (Case 4). Here, the (7 -) 

K sr2 max 


value is 1.80, which is 4% 


greater than that for the unstiffened one. 
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IV. Concurrent Processing 


All of the global analysis computations for the aluminum and composite panels 
investigated in this study were carried out on a parallel computer (Flex/32 Multicomputer, 
reference 2.) installed at NASA Langley Research Center. For each problem, the global 
system of equations for finding the nodal displacements was solved on up to 16 parallel 
processors using the Choleski method (reference 1). Table 15 gives the computational times 
and speedups for the unstiffened and stiffened aluminum panels with hole. Table 16 gives 
similar results for the aluminum panel with hole and inner or outer stiffeners. The 
computational time and speedup results for the unstiffened and stiffened panels without hole 
was found to be quite similar to that given in Table 15 for the panels with hole. Similarly, 
the results for the composite panels were found to be practically the same as that for the 
aluminum panels. As illustrations. Figures 21 and 22 present the relationships between the 
computational time (T c ) versus the number of parallel processors corresponding to the results 
given in Table 15. The associated speedup curves are plotted in Figure 23 and compared to 
the ideal relationship. The advantage of concurrent processing for solving finite element 
analysis problems of the type considered herein is quite obvious from these figures. 

The bandwidths for the global stiffness matrix for the panels shown in Figures 1(c) 
through 1(f) are found to be 70, 195, 218, and 240, respectively. Figure 24 shows the 
variation of the computational time T c , given in Tables 15 and 16, with matrix bandwidth B m , 
and the number of processors N. The curves in this figure indicate that the steepness of the 
T c -B m relationship decreases with an increase in N. This may be attributed to increased 
inter- processor communication time losses. 
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V. CONCLUSIONS 


The study of the distribution of stresses for unstiffened and stiffened, aluminum and 

composite panels has led to the following useful and interesting conclusions: 

1. The maximum compressive stresses for composite panels are substantially smaller than 
those for aluminum panels. 

2. Local analysis of panels provides a detailed and accurate distribution of stresses 

concentrated around the hole. 

3. The outer stiffeners in the aluminum panel develop higher compressive stresses than 
does the inner one. This difference is even greater in the composite panel. 

4. In both aluminum and composite panels, the maximum shear stress exhibits a greater 
variation along the panel length at the inner stiffener location in comparison to that 
observed for the outer stiffener location. 

5. The maximum compressive stresses from the global analysis for both stiffened aluminum 

and composite panels are greater than the corresponding stresses for the unstiffened 

ones. The panels with only the inner stiffener develop greater stresses than those for 
the unstiffened panels. However, the outer two stiffeners cause a reduction in the stress 
values, although this reduction does not completely off-set the detrimental influence of 
the inner stiffener. 

6. The layer-by-layer stress analysis shows that the maximum value of the principal 

compressive stress exists on the 0* ply-orientation for both stiffened and unstiffened 
composite panels. The maximum principal tensile stress exists on the -45° ply-orientation 
for the unstiffened panel, and on the bottom-most 45° ply-orientation for the stiffened 
panel. The largest value of the maximum shear stress occurs on the 0° ply-orientation 
for the unstiffened composite panel, and on the bottom most 45' ply-orientaion for the 
stiffened composite panel. 

7. The computational time required by the Choleski method decreases as the number of 
processors is increased and decreases with the bandwidth. 
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Table 1. Boundary Conditions for Unstiffened Panels 
(u=unconstrained; c=constrained) 
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Table 2. Boundary Conditions for Stiffened Panels 
(u=unconstrained; c=constrained) 
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Stiffeners 

Deflection 

Edges 

Corners 

Edges 


OP 

PQ 

QS 

SO 

0 

D 

Q 

s 

KO IJ 

AXj 

u 

c 

u 

u 

u 

c 

c 

u 

u c 

AX2 

u 

u 

u 

u 

u 

u 

c 

c 

c c 

ax 3 

u 

c 

u 

c 

c 

c 

c 

c 

c c 

*Xj 

u 

c 

u 

c 

c 

c 

c 

c 

c c 

« 2 

u 

c 

u 

c 

c 

c 

c 
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c c 

*X 3 
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c 

u 

c 

c 

c 

c 

c 

c c 
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Table 3. Dimensionless Principal Stresses, 02 , for Aluminum 
Panels without Hole. 


Location 

(inches) 

Dimensionless Principal Stress, 02 



Unstiffened 

Stiffened 

*1 

x 2 





(33333.2) 

(33918.8) 

0.00 

0.00 

1.00 

1.00 

0.00 

5.75 

1.00 

0.99 

0.00 

11.50 

1.00 

1.02 

15.00 

0.00 

1.00 

0.98 

15.00 

5.75 

1.00 

0.98 

15.00 

11.50 

1.00 

0.98 


* o 0 » P*» 
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Table 4. Dimensionless Principal Stresses, * 2 , for Stiffeners 
in Almpinum Panels without Hole 


Location 

(inches) 

Dimensionless Principal 
Stress, *2 

Xj 

x 2 

x 3 

(33918.80)* 

0.00 

1.25 

0.00 

0.99 

0.00 

1.25 

1.40 

0.99 

15.00 

1.25 

0.00 

0.99 

15.00 

1.25 

1.40 

0.99 

0.00 

5.75 

0.00 

0.99 

0.00 

5.75 

1.40 

0.99 

15.00 

5.75 

0.00 

0.98 

15.00 

5.75 

1.40 

0.99 

0.00 

10.50 

0.00 

0.99 

0.00 

10.50 

1.40 

0.98 

15.00 

10.50 

0.00 

0.98 

15.00 

10.50 

1.40 

0.99 


ffQ.PSl 
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Table 5. Dimensionless Principal Stress Resultants, 5 ^, for 
Composite Panels without Hole. 


Location 

(inches) 

Dimensionless Principal Stress 
Resultant, 5^ 



Unstiffened 

Stiffened 

Xl 

*2 

(5135.87)* 

(5344.77) 

0.00 

0.00 

1.00 

1.00 

0.00 

5.75 

1.00 

0.96 

0.00 

11.50 

1.00 

1.05 

15.00 

0.00 

1.04 

1.00 

15.00 

5.75 

1.13 

1.05 

15.00 

11.50 

1.04 

1.01 









Table 6. Dimensionless Principal Stress Resultants, ?sr2* for 
Stiffeners in Composite Panels without Hole ^ 


Location 

(inches) 

Dimensionless Principal 
Stress Resultant 

Xl 

*2 

*3 

(5344.77)* 

0.00 

1.25 

0.00 

1.10 

0.00 

1.25 

1.40 

1.39 

15.00 

1.25 

0.00 

1.10 

15.00 

1.25 

1.40 

1.33 

0.00 

5.75 

0.00 

1.10 

0.00 

5.75 

1.40 

1.38 

15.00 

5.75 

0.00 

1.12 

15.00 

5.75 

1.40 

1.20 

0.00 

10.50 

0.00 

1.09 

0.00 

10.50 

1.40 

1.33 

15.00 

10.50 

0.00 

1.17 

15.00 

10.50 

1.40 

1.37 









Tabic 7. Maximum principal stresses around hole for various 
layers of unstiffened composite panel (psi) 


5055?. 37 
45 50553.37 

0 1 1 1 1 5.97 

0 11115.97 

-45 50553.37 

45 5055?. 37 

0 11115.97 

0 11115.97 

0 11115.97 

45 5055?. 37 

-45 50553.37 

0 11115.97 

11115.97 
11115.97 
50553.37 
5055?. 37 
11115.97 
11115.97 
11115.97 
5055?. 37 
50553.37 
11115.97 
11115.97 
50553. 37 
5055?. 37 


2 

3 

4 

5 

6 

7 

8 

.9 

10 

II 

32970.15 

17456.16 

11925.35 

-1621.11 

-132 3.72 

333.71 

19960.84 

58137.94 

3402 7.39 

20594.74 

38267.97 

1423.72 

-392.66 

-1613.02 

10? 19. 38 

20594. 35 

34026.01 

58139. 24 

19963 . 06 

338. 88 

8028.16 

1816.08 

196.70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

640 7.87 

-268.88 

8028. 16 

1816.08 

196. 70 

914.88 

-1116. 95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

38267.97 

1423.72 

-392.66 

-1613.02 

10219.38 

20594.35 

34026.01 

58139.24 

19963.06 

338.88 

32970.15 

17456.16 

11925.35 

-1621.11 

-1323. 72 

333. 71 

19960.84 

58137.94 

34027.39 

20594. 74 

8028.16 

1816.08 

196.70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

8028. 16 

1816.08 

1%. 70 

914.88 

-1116. 95 

-269.23 

6407.55 

12488.00 

6407.87 

-268. 88 

8028.16 

1816.08 

196.70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

32970.15 

17456.16 

11925.35 

-1621.11 

-1323. 72 

333. 71 

19960.84 

58137.94 

34027.-39 

20594. 74 

38267.97 

1423.72 

-392.66 

-1613.02 

10219.38 

20594.35 

34026.01 

58139.24 

19963.06 

338.88 

8028.16 

1816.08 

1%. 70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268. 88 

8028.16 

1816.08 

196.70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

8028.16 

1816.08 

1%. 70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

38267.97 

1423.72 

-392.66 

-1613.02 

10219.38 

20594.35 

34026.01 

58139.24 

19963.06 

338.88 

32970.15 

17456.16 

11925.35 

-1621 . It 

-1323. 72 

333. 71 

19960. 84 

58137.94 

34027. 39 

20594. 74 

8028.16 

1816.08 

196.70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

8028.16 

1816.08 

1%. 70 

914.88 

-1116. 95 

-269.23 

6407.55 

12488. 00 

6407.87 

-268. 88 

8028.16 

1816.08 

196.70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

32970. 15 

1 7456. 16 

11925.35 

-1621.11 

-1323. 72 

333. 71 

19960.84 

58137.94 

34027.39 

20594. 74 

38267.97 

1423.72 

-392.66 

-1613.02 

10219.38 

20594.35 

34026.01 

58139.24 

19963.06 

338.88 

8028. 16 

1816.08 

1 96. 70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

8028.16 

1816.08 

196.70 

914.88 

-1116.95 

-269.23 

6407.55 

12488.00 

6407.87 

-268.88 

38267.9? 

1423. 7? 

-392.66 

-1613.02 

10219.38 

20594.35 

34026.01 

58139.24 

19963.06 

338.88 

32970.15 

17456.16 

11925.35 

-1621.11 

-1323.72 

333.71 

19960.84 

58137.94 

3402 7.39 

20594.74 


10 ? 1 7.73 
-1 3?3. 75 
-1116.77 
-1116.77 
-13? 3.75 
10?1 7. 73 
-1116.77 
-1116.77 
-1116.77 
10?1 7. 73 
-1323.75 
-1116.77 
-1116.77 
-1116.77 
-1323.75 
10217.73 
-1116.77 
-1116.77 
-1116.77 
10217.73 
-1323.75 
-1116.77 
-1116.77 
-1323. 75 
10217.73 


-1613.16 -392.66 1423.59 38265.45 45 

-1621.19 11926.19 17457.15 32969.35 -45 

914.76 196.63 1815.86 8027.86 0 

914.76 196.63 1815.86 8027.86 0 

-1621.19 11926.19 17457.15 32969.35 -45 

-1613.16 -392.66 1423.59 38265.45 45 

914.76 196.63 1815.86 8027.86 0 

914.76 196.63 1815.86 8027.86 0 

914.76 196.63 1815.86 802 7.86 0 

-1613.16 -392.66 1423.59 38265.45 45 

-1621.19 11926.19 17457.15 32969.35 -45 

914.76 196.63 1815.86 8027.86 0 

914.76 196.63 1815.86 8027.86 

914.76 196.63 1815.86 8027.86 

-1621.19 11926.19 17457.15 32969.35 

-1613.16 -392.66 1423.59 38265.45 

914.76 196.63 181 S .86 8027.86 

914.76 196.63 1815.86 8027.86 

914.76 196.63 1815.86 8027.86 

-1613.16 -392.66 1423.59 38265.45 

-1621.19 11926.19 17457.15 32969.35 

914.76 196.63 1815.86 8027.86 

914.76 196.63 1815.86 8027.86 

-1621.19 11926.19 17457.15 32969.35 

-161 3.16 -392.66 1423.59 38265.45 
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Table 8. Minimum principal stresses around hole for various 
layers of unstiffened composite panel (psi) 


r 

jo ^ 

>| 


k 

i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

.■-< 0q 

14 15 

16 

J / 

i 

5364.10 

3736.57 

- 1761.31 

- 6094.23 

-29115.46 

-73510.96 

-52004.36 

2582.18 

6345.70 

2280.81 

- 4613.85 

- 8091.88 

-28478.57 

-53938.37 

-19656.38 

3556.98 

45 

-45 

5364.05 

3557.17 

-19652.84 

-53936.12 

-28476.72 

- 8092.01 

- 4614.09 

2280. 54 

6345.65 

2582.49 

-52001.54 

-73506.98 

-29113. 79 

- 6094.35 

- 1761.57 

3736. 30 

-45 

0 

-2213/.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049.75 

-39102.79 

-7/686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

0 

-22127.42 

-18966.05 

-30253. 34 

-66397.49 

-96564. 11 

-77686.97 

-39101.51 

-20050. 32 

-20352. 12 

-20049. 75 

-39102. 79 

-77686. 88 

-96562.81 

-66398.24 

-30253.42 

-18966.6 1 

0 

-45 

5364.10 

3557.17 

-19652.84 

-53936.12 

-28476.72 

- 8092.01 

- 4614.09 

2280.54 

6345.65 

2582.49 

-52001.54 

-73506.98 

-29113.79 

- 6094.35 

- 1761.5/ 

3736.30 

-45 

45 

5364. 10 

3736.57 

- 1761.31 

- 6094.23 

-29115.46 

-73510.96 

-52004. 36 

2582. 18 

6345.70 

2280.81 

- 4613.85 

. - 8091.88 

-28478.57 

-53938.37 

-19656.38 

3556.98 

45 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049.75 

-39102.79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

0 

-22127.42 

-18966.05 

-30253. 34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050. 32 

-20352.12 

-20049. 75 

-39102. 79 

-77686. 88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049.75 

-39102.79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

u 

45 

5364.10 

3736.57 

- 1761.31 

- 6094.23 

-29115.46 

-73510. 96 

-52004.36 

2582. 18 

6345. 70 

2280.81 

- 4613.85 

- 8091.88 

-28478.57 

-53938.37 

-19656.38 

3556.98 

45 

-45 

5364.10 

3557.17 

-19652.84 

-53936.12 

-28476.72 

- 8092.01 

- 4614.09 

2280.54 

6345.65 

2582.49 

-52001.54 

-73506.98 

-29113.79 

- 6094.35 

- 1761.57 

3736.30 

-45 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049. 75 

-39102. 79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

t> 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049.75 

-39102.79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050. 32 

-20352.12 

-20049. 75 

-39102. 79 

-7 7686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

-45 

5364.10 

3557.17 

-19652.84 

-53936.12 

-28476.72 

- 8092.01 

- 4614.09 

2280.54 

6345.65 

2582.49 

-52001.54 

-73506.98 

-29113.79 

- 6094.35 

- 1761.57 

3736.30 

-45 

45 ' 

5364.10 

3736.57 

- 1761.31 

- 6094.23 

-29115.46 

-73510.96 

-52004.36 

2582. 18 

6345. 70 

2280.81 

- 4613.85 

- 8091.88 

-28478.57 

-53938.37 

-19656.38 

3556.98 

45 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049.75 

-39102.79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

0 

-22127.42 

-18966.05 

-30253. 34 

-66397.49 

-96564. 11 

-77686.97 

-39101.51 

-20050. 32 

-20352. 12 

-20049.75 

-39102. 79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049.75 

-39102.79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

45 

S364. 10 

3736.57 

- 1761.31 

- 6094.23 

-29115.46 

-73510.96 

-52004. 36 

2582. 18 

6345. 70 

2280. 81 

- 4613.85 

- 8091.88 

-28478.57 

-53938.37 

-19656.38 

3556.98 

45 

-45 

5364.10 

3557.17 

-19652.84 

-53936.12 

-28476.72 

- 8092.01 

- 4614.09 

2280.54 

6345.65 

2582.49 

-52001.54 

-73506.98 

-29113.79 

- 6094.35 

- 1761.57 

3736.30 

-4$ 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564. 11 

-77686.97 

-39101.51 

-20050. 32 

-20352. 12 

-20049. 75 

-39102. 79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

0 

0 

-22127.42 

-18966.05 

-30253.34 

-66397.49 

-96564.11 

-77686.97 

-39101.51 

-20050.32 

-20352.12 

-20049.75 

-39102.79 

-77686.88 

-96562.81 

-66398.24 

-30253.42 

-18966.63 

(J 

-45 

5364. 10 

3557.1 7 

-19652.84 

-53936.12 

-28476. 72 

- 8092.01 

- 4614.09 

2280. 54 

6345.65 

2582.49 

-52001.54 

-73506.98 

-29113. 79 

- 6094.35 

- 1761.57 

3736. 30 

-45 

45 

5364.10 

3736.57 

- 1761.31 

- 6094.23 

-29115.46 

-73510.96 

-52004.36 

2582.18 

6345.70 

2280.81 

- 4613.85 

- 8091.88 

-28478.57 

-53938.37 

-19656.38 

3556.98 

4:* 

n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

II 

12 

13 

14 

15 

16 

\tf 

j\ 









Table 9. Maximum shear stresses around hole for various layers 
unstiffened composite panel (psi) 


X 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

X 



45 

22594.13 

14616.79 

9608.73 

9009.79 

1374 7.17 

36093.62 

26169.03 

8689.33 

25896.1 2 

15873.29 

12604.30 

9154.80 

13432.71 

26772.85 

10539.98 

17354.24 

45 



-45 

22594.66 

17355.40 

10538.28 

26771. 73 

13431.85 

9155. 70 

12604.22 

15872. 73 

25896. 80 

8690.28 

26167. 71 

36091.61 

13746.30 

9010.2 7 

9609. 36 

14616.52 

-45 



0 

16621.70 

13497.10 

16034.71 

33297.09 

48739.50 

38285.01 

19416.14 

13228.93 

16420.06 

13228.81 

19416.96 

38285.05 

48738.79 

33297.44 

16034.64 

13497.25 

0 



0 

16621. 70 

13497.10 

16034. 71 

33297.09 

48739.50 

38285.01 

19416.14 

13228. 93 

16420. 06 

13228.81 

19416.96 

38285. 05 

48738.79 

33297.44 

16034.64 

1 3497. 25 

0 



-45 

22594.70 

17355. to 

10538.28 

26771.73 

13431.85 

9155.70 

12604.22 

15872.73 

25896.80 

8690.28 

26167.71 

36091.61 

13746.30 

9010.27 

9609.36 

14616.52 

-45 



45 

22594.10 

14616.79 

9608. 73 

9009. 79 

13747.17 

36093.62 

26169.03 

8689.33 

25896. 12 

15873.29 

12604. 30 

9154.80 

13432. 71 

26772.85 

10539.98 

17354.24 

45 



0 

16621.70 

13497.10 

16034.71 

33297.09 

48739.50 

38285.01 

19416.14 

13228.93 

16420.06 

13228.81 

19416.96 

38285.05 

48738.79 

33297.44 

16034.64 

13497.25 

0 

O 

O 




















5*3 

0 

16621. 70 

13497.10 

160 34 . 71 

33297.09 

48739.50 

38285.01 

19416.14 

13228.93 

16420. 06 

13228.81 

19416.96 

38285. 05 

48738. 79 

33297.44 

16034.64 

13497.25 

0 

►d 

O 

0 

16621.70 

13497.10 

16034.71 

33297.09 

48739.50 

38285.01 

19416.14 

13228.93 

16420.06 

13228.81 

19416.96 

38285.05 

48738.79 

33297.44 

16034.64 

13497.25 

0 

o 

o 

y 

45 

22594.10 

14616.79 

9608. 73 

9009. 79 

13747.17 

36093. 62 

26169.03 

8689. 33 

25896.12 

15873.29 

12604. 30 

9154.80 

13432. 71 

26772.85 

10539. 98 

17354.24 

45 

5*3 

> 

r 1 , 

-45 

22594.70 

17355.40 
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Table 10. Maximum principal stresses around hole for various 
layers of stiffened composite panel (psi) 
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Table 11 . Minimum principal stresses around hole for various 
layers of stiffened composite panel (psi) 
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Tabic 12. Maximum shear stresses around hole for various layers 
stiffened composite panel (psi) 
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Table 13. Effect of Stiffeners on Stress Concentration Factor for 
Aluminum Panels with hole based on Global Analysis 



Table 14. Effect of Stiffeners on Dimensionless Maximum Stress Resultants 
for Composite Panels with hole based on Global Analysis 
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Table 15. Computational Times and Speedups for Unstiffened and Stiffened 
Aluminum Panels with Hole 


Number 

of 

Processors 

Computational Time 
(seconds) 

Speedup " 

Unstiffened 

Stiffened 

Unstiffened 

Stiffened 

1 

22.2 

684.0 

_ 

- 

2 

11.4 

346.0 

1.9 

2.0 

4 

6.2 

179.0 

3.6 

3.8 

8 

4.3 

92.3 

5.1 

7.4 

16 

4.1 

50.5 

5.4 

13.5 


Table 16. Computational Times and speedups for Aluminum Panel with Hole 
and Inner Stiffener (I.S.) or Outer Stiffener(O.S.) 


Number 

of 

Processors 

Computational Time 
(seconds) 

Speedup 

with I.S. 

with O.S. 

with I.S. 

with O.S. 

1 

408.0 

486.0 

- 

- 

2 

206.0 

245.0 

2.0 

2.0 

4 

106.0 

127.0 

3.8 

3.8 

8 

55.7 

65.7 

7.3 

7.4 

16 

31.6 

36.6 

12.96 

13.3 
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(c) Unstiffened Panel with Hole 


(f) Stiffened Panel with Hole 


Figure 1. Types of panels 
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(b) Stlfftatr Dlieretlaatloa 



(d) Global/Local Diicretlialion Arouod Hole 



(«) Dltcretiiatfoa for Stiffened Panel with Hole 


Figure 2. Finite Element Discretizations 
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(b) Outer stiffener 



Figure 4. Contours for dimensionless compressive principal stress for stiffened aluminum panel 
by global analysis (<T 0 _ 32994.29 psi) 
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Figure 5. Midsurface maximum shear stress distribution for aluminum panel at location of 
stiffeners 
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Figure 8. Distribution of maximum shear stress resultants for composite panel at location of 
stiffeners 
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rincipal stresses for various layers of stiffened composite panel 
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Figure IS. Minimum principal stresses for various layers of 
stiffened composite panel 
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Figure 23. Speedup curves for unstiffened and stiffened aluminum 
panels with hole 
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